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We use the (5iV-formalism to describe the leading order contributions to the primordial power 
spectrum, bispectrum and trispectrum in multiple-field models of inflation at leading order in a 
, perturbative expansion. In slow-roll models where the initial field fluctuations at Hubble-exit are 

■ nearly Gaussian, any detectable non-Gaussianity is expected to come from super-Hubble evolution. 

' We show that the contribution to the primordial trispectrum can be described by two non-linearity 

£\J , parameters, tjvl and gNL, which are dependent upon the second and third derivatives of the local 

^ . expansion with respect to the field values during inflation. 

q| PACS numbers: 98.80.Cq PU-ICG-06, |astro-ph/0611075| 
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The simplest models of inflation predict a quasi scale invariant spectrum of nearly Gaussian, adiabatic perturba- 
tions While observations are consistent with this picture, it is important to measure any deviations from these 
£ — \ predictions, in order to constrain the many models of inflation. Non-Gaussianity is a potentially powerful discriminant 
■ between different models. Currently most efforts to constrain the non-Gaussianity have focused on the 3-point func- 
tion of perturbations, the bispectrum 0J3Jj|. However the 4-point function, or trispectrum can also be constrained 
by increasingly accurate measurements [5|, |f| |7j . 

The JiV-formalism [H, [^, ITol . [Tl| identifies the primordial curvature perturbation with a perturbation in the local 
value of the integrated expansion, 
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O 

ryj , with respect to the expansion in the background spacetimc 



N(x) = I " H(x,t)dt, (1) 



N = I " H(t)dt . (2) 



^ . This enables one to calculate the non-linear primordial curvature perturbation (e.g., at the epoch of primordial 
nucleosynthesis) in terms of initial scalar field fluctuations during inflation [ll[ ■ In particular the primordial curvature 
perturbation on uniform-density hypcrsurfaccs, £, corresponds to the perturbation in the local expansion defined 
with respect to an initial spatially flat hypersurface [To| . In the long- wavelength limit where spatial gradients and 
anisotropic shear becomes small, the integrated expansion along a worldlinc can be calculated from solutions to the 
unperturbed Fricdmann equation. 
We have 

( = N(0 A )-N, (3) 

where (f> A represents the initial values for the scalar fields on the initial spatially-flat hypersurface. If we decompose 
the local scalar field values into a homogeneous background value plus perturbation, 

4>A = <t>A + VA , (4) 
we can write the curvature perturbation as a Taylor expansion 

(=E%-^V 2 .../-. (5) 
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There is an implicit sum over the fields for all of the repeated indices A^. We are thus able to describe the non- 
Gaussianity of the primordial curvature perturbation given an initial distribution of scalar field perturbations. In 
practice, scalar fields during slow-roll inflation arc expected to have a nearly Gaussian distribution during inflation, at 
least on scales slightly larger than the Hubble scale, i.e., shortly after Hubble exit [l2||. In single-field inflation this is 
sufficient to guarantee that the non-Gaussianity of the primordial curvature perturbation is very small (l3l . Q, [H| ■ It is 
suppressed by slow-roll parameters since one can express derivatives of the expansion in terms of slow-roll parameters. 
However in multi-field models the perturbed expansion cannot necessarily be related to slow-roll parameters at Hubble- 
exit, and the primordial curvature perturbations may deviate strongly from a Gaussian distribution. 



II. PERTURB ATIVE EXPANSION 



In linear perturbation theory the initial free-field fluctuations describe Gaussian random fields, which we will denote 
by ^>i{t, x % ). In a perturbative expansion higher-order interactions and non-linear evolution lead to higher-order terms 
which are non-Gaussian 

<f A = ft + \v>2 + + ■ ■ ■ • (6) 

From the definition of Gaussian statistics, the n-point correlator of the Gaussian fluctuation is zero for odd n and 
can be reduced to (disconnected) products of two-point functions for even n. 



A. Power spectrum 



In Fourier space the power spectrum is given by 

(^k'>=^ B (fc)(2^ 3 (k + k'). 



(7) 



This is second and higher order in a perturbative expansion. 

At leading order in the field perturbations and in the slow-roll limit the fluctuations are independent and we have 

C AB {k) = S AB P(k), (8) 

where 6 AB is the Kronecker delta-function, and the variance per logarithmic interval in fc-space is given by 



V(k) = 



4irk 3 



I If 



(9) 



where the Hubble parameter H is evaluated at Hubble-cxit, k = (ai?)*. At zcroth order in slow-roll parameters, V is 
independent of wavenumber, i.e. we have a scale invariant spectrum for the field fluctuations. 

In general the fields are correlated at Hubble-exit, at first order in the slow-roll parameters. In the case of two field 
inflation, using the methods of (see also we find for the power spectra and cross-correlation respectively, 
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where C = 2 — In 2 — 7~ 0.7296 and 7 is the Eulcr-Maschcroni constant. Overdots represent derivatives with respect 
to time, and V < f >1 = dV/d(f>i. 
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B. Bispectrum 

The first signal of non-Gaussianity comes from the bispectrum, which at lowest order is 

(<p A <p B <p C ) = {<p£<p?<p?) + perms , (13) 

and therefore is fourth order in perturbations. 
The bispectrum of the distribution is given by 

« < <£) = B ABC (k u k 2 , fc 3 )(2^) 3 ( 5 3 (k 1 + k 2 + ka) . (14) 

This is fourth and hi ghe r order. This was originally calculated by Maldacena, [l3[ for single field inflation, and 
by Seery and Lidsey, [l2|, for multiple fields. They show that it only depends on the amplitude of the k vectors. 
Specifically they calculate the quantity A ABC at leading order in slow roll, which is related to B ABC by 

4-TT 4 

B ABC {k u k 2 , h) = ^V 2 A ABC {k u k 2 , k 3 ) . (15) 
They find A ABC ~ 0(e 1 ^ 2 ) so it vanishes in the slow-roll limit. 



C. Trispectrum 

From Wick's theorem, the first-order, Gaussian perturbations do not contribute to the connected part of the four- 
point function, 

(tfyV?pf>e = 0. (16) 

Note there is a disconnected part of the four-point function which is a product of two two-point functions which is 
only fourth order in slow roll and arises even for purely Gaussian statistics. This disconnected term is only non-zero 
if, e.g., the momenta satisfy ki + k2 = and k 3 + k 4 = 0. 

The second-order field perturbations, ip A , are generated from the product (or the convolution in Fourier space) of 
two first order variables that have Gaussian distributions, so they do not contribute to the (connected) four-point 
function at lowest possible order, 

fopMV?P?>c = 0. (17) 
Hence the leading order contribution to the connected four-point function comes from two terms, [f|, 

(<Pi<Pi<Pi<P3) c and (vifi^^c . (18) 

Both of these terms arc sixth-order in perturbations. 

In Fourier space the connected part of the four-point function is sixth order, and given by 

(<<<^>k 4 )c = T' 4SCI? (k 1 ,k 2 ,k3,k 4 )(2^ 3 ( 5 3 (k 1 +k 2 +k3+k 4 ). (19) 

This was recently calculated at leading order by Seery, Lidsey and Sloth [l8[. They show this quantity depends on 
both the magnitude and direction of the k vectors, specifically it depends on ki = |kj| and k^ ■ kj. 



III. THE PRIMORDIAL N-POINT FUNCTIONS 



So far we have calculated the two-, three- and four-point function of the field fluctuations. To link these to 
observations, we need to calculate the n-point functions of the primordial curvature perturbation £. We do this using 
the 6N expansion for (, ([5|). 
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A. The primordial power spectrum 

At leading order the primordial power spectra depends purely on £i, from Eq. (O, 

Ci = N A cp A . (20) 

The power spectrum is thus 

(CkCk<) =P c (k)(27rfS 3 (k + k'), (21) 

where 

P c (k) = N A N B C AB (k). (22) 

In the slow-roll limit, ((8|), this reduces to 

P c (k) = N A N A P(k) . (23) 

B. The primordial bispectrum 

To leading order in the field perturbations, the 3-point function of the curvature perturbations depends on £i, (|20p 
and 

C2 = N A ip A + N AB ip A <pf . (24) 

The primordial bispectrum is thus 

(C kl Ck 2 Ck 3 ) = N A N B N G (^i <p B 2 ) 

+\n MM N b N c [(& Al * V M )^ + (2 perms)] , (25) 

where V denotes the convolution, defined by 

(^*^)k = ^/^*V£- k '*£. (26) 

Hence the bispectrum of the curvature perturbation is 

(Ck, Ck 2 Ck 3 ) = B c {k u ka, fc 3 )(27r) 3 5 3 (k 1 + k 2 + k 3 ) , (27) 

where to leading order [l9| 

B c (k u k 2 ,k 3 ) = N A N B NcB ABC (ki,k 2 ,k 3 ) 

+N A N BC N D [C AC (h)C BD (k 2 ) + C AC {k 2 )C BD {h) + C AC {k z )C BD {k 1 )} . (28) 



In the slow-roll limit we can write the bispectrum as [2C 



A^kf^f -1 HkiMM) N AB N A N l 



where the form factor T is defined by 



{NcNcf I ' 



(29) 



H ku k 2 ,k 3 ) = £ kf - £ hk) - S ^^l • (30) 
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C. The primordial trispectrum 

From the discussion in sec. Ill CI the four-point function of the curvature perturbation at leading order will depend 
on Ci, ©, Ca, HMD, and 

Cs = N A <fli + N AB (<p A ip B + <p A ip B ) + N ABC <p A <p?ip? (31) 
The four-point function at leading order is 

(Ck 1 Ck 2 Ck 3 Ck 4 ) c = N A N B N c N D (ip^<p^ 3 ip^) c 

+ ]-N AlA2 N B N c N D \((<p A > * ^ 2 ) ki + (3 perms) 



2 
1 

T 

, ^N AlA2A3 N B N c N D [((^ * ip A * * ^ 3 ) kl + (3 Perms)] . (32) 

All of the terms shown are sixth order in the field perturbations. 

The first term of the expansion above is the intrinsic 4-point function of the fields, as calculated, [3]. The 
disconnected part of this term would only give a contribution if the sum of any two k vectors is zero, e.g. if ki +k2 = 0. 
We will exclude this case, which is equivalent to neglecting parallelograms of the wavevectors. 

The second term of (|3"!?)) consists of permutations of terms of the form 

1 -N MA2 N B N C N D ^ J A«< 2 _XX 3 <) (33) 

This five-point function is zero for the first-order, Gaussian, perturbations, hence the leading order contribution is 
sixth-order, due to the second order contribution of one of the fields. Hence we use Wick's theorem to split the 5-point 
function in to lower point functions. There is no contribution to (|33p from the split into a four-point and a one-point 
function. Only the split into a two-point and three-point function gives a contribution. However the first possible 
contraction in (|33|) . (Phi-a^t 2 ) > does not contribute since it is only non-zero when k\ = 0. Therefore we can reduce 
the above term into a power spectra and a trispectrum in 6 different ways, which gives three distinct pairs of terms. 
In total the second term of i 



+ -N AlA2 N BlB2 N c N D <(^ Al * v^ 2 ) kl (<P Bl * ^ S2 ) k2 + (5 perms) 



N AlA2 N B N c N D [C AlS (fc 1 )B A2BC (fc 12 ,fc 3 ,fc 4 ) + (11 perms)] (2^) 3 <5 3 (k t ) , (34) 

where we use the shortened notation kij = |k» + kj| and kt = ki + k 2 + k 3 + k 4 . The 12 permutations come from 
having 3 distinct choices for the indices of the wavenumber kij (only three distinct choices because kij = kji and 
ki2 = &34 etc). We then choose which two wavenumbers form the remaining arguments of B ABC , either ki, kj or the 
other pair of wavenumbers, and finally we choose which of the two indices i or j is attached to the wavenumber ki 
that is the argument of C. 

Continuing this argument for the second and third terms of (|32[) . we find the connected part of the trispectrum of 
the curvature perturbation is 



(Cka Ck 2 Ck 3 Ck 4 )c = T c (ki, k 2 , k 3 , k 4 )(2 7 r) 3 ( 5 3 (k 1 + k 2 + k 3 + k 4 ) , (35) 



where 



^(ki.ka.ka.k*) = A^A^cA^T^^k^^k^k 



+N AlA2 N B N c N D [C MB {ki)B^ BC {kvi,kiM + (11 perms)] 
+N AlA2 N BlB2 N c N D [C A2B2 (k 13 )C AlC (k 3 )C BlD (h) + (11 perms)] 
+N AlA2A3 N B NcN D [C MB (k 2 )C MC {k 3 )C MD {k i ) + (3 perms)] . (36) 

IV. GAUSSIAN SCALAR FIELDS 

If the scalar field perturbations are independent, Gaussian random fields, as we expect shortly after Hubble-cxit 
during inflation in the slow-roll limit [HI, EH then the bispectrum for the fields, B ABC , and connected part of the 
trispectrum, T ABCD , both vanish. 
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/^= r,.: 2 ■ (38) 



In this case the bispcctrum of the primordial curvature perturbation (|28[) at leading (fourth) order, can be written 

as 

S f (kx,k 2 ,k 3 ) = lfNL[Pdk 1 )Pc(k2)+Pdk2)Pdh)+P c (k 3 )P c (k 1 )} . (37) 

where the dimensionless non-linearity parameter 1 is given by [TlT ] 

5 N A N B N AB 

6 (N c N°f 
The trispectrum ([36]) in this case reduces to 

T c (k 1; k 2 ,k3,k 4 ) = N AB N AC N B N C [P(k 13 )P(k 3 )P(k 4 ) + (11 perms)] 

+N ABC N A N B N C [P(k 2 )P(k 3 )P(k 4 ) + (3 perms)] , (39) 

Hence we can write the trispectrum as 

T c (ki,k 2 ,k3,k4) = t nl [P c (fc 13 )P c (fc 3 )P c (fc4) + (11 perms)] 
54 

+ ^.9jvl [P C ( W(W( fc 4) + (3 perms)] . (40) 
where comparing the above two expressions, and using (|23[) we see 

N AB N AC N B N C 

™ L = (N D N°Y ' (41) 

57Vi = 54 (iV c iV^)3 • (42) 

The expression for tjvl from multiple fields was given in the arXiv version of (2lj . Note that we have factored 
out products in the trispectrum with different k dependence in order to define the two k independent non-linearity 
parameters tnl and qnl- This gives the possibility that observations may be able to distinguish between the two 
parameters 



A. Single field dependence 

In many cases there is single direction in field-space, <p, which is responsible for perturbing the local expansion, 
N(ip), and hence generating the primordial curvature perturbation, £. For example this would be the inflaton field in 
single field models of inflation, or it could be the late-decaying scalar field in the curvaton scenario [22j . 

In the case where a single field dominates, the curvature perturbation ([5]) is given by 

C = N'ip + \N"^ + VV + • • • , (43) 

where we use the shorthand N' — dN/dip. If in addition we assume that the field perturbation is purely Gaussian, 
(p — ipi, then the non-Gaussianity of the primordial perturbation has a simple "local form" where the full non-linear 
perturbation at any point in real space, C( x )i is a local function of a single Gaussian random field, ip±. Thus we can 
write 0, [H 

C = Ci + IfNLC! + ^9nlCi +■■■, (44) 

where Ci is Gaussian because it is directly proportional to the initial Gaussian field perturbation, ipi, and the dimen- 
sionless non-linearity parameters, //vz, and g^L, ar e given by 

5 N " fAKS 

fNL = 6(iV^' ( } 

25 N "' fABS 

9NL = 54 (W ( } 



1 Some papers use a different sign convention for fjfL- F° r example, Refs. lilt |20H use the opposite sign convention. 
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The numerical factors in Eq. (|44[) arise because the original definition is given in terms of the Bardeen potential on 
large scales (in the matter dominated era, md), &Hmd = (3/5)Ci, so we have 0,0,0], 

3 

-C = ®Hmd + fNL&HmA + 9NL$Hmd + ' ' ' ■ (47) 



The primordial bispcctrum and trispectrum are then given by Eqs. (|3T[) and (|40|) . where the non-linearity parameters 
Jnl and gNL, given in Eqs. (|38|) and (|42|) . reduce to Eqs. (|45| and (146]) respectively, and tatl given in Eq. (j4l"j) reduces 
to 

NL = Jn 7 ) 1 ~zh iNL ' ( ' 

Notice that the bispcctrum depends linearly on Q2 while the trispectrum has both a quadratic dependence upon £2 
and a linear dependence on £3. Thus tnl is proportional to f NL (shown in 0] using the Bardeen potential, and in 
[24| using this notation). However the trispectrum could be large even when the bispectrum is small because of the 
gNL term 0, [23|. 



1. Inflaton scenario 



In the case of standard single held inflation, where the primordial curvature perturbation is generated solely by 
the inflaton field, we can calculate the non-linearity parameters f^L and g^L in terms of the slow-roll parameters at 
Hubblc-cxit. Because the large-scale perturbations are adiabatic, C is non-linearly conserved on large scales [Tol.[25l.[2rl] 
and the derivatives of the expansion, N' , N" and N'" can be calculated at Hubble-exit. Using the definition (JTJ) , we 
find 

N> = L^I^OU) , (49) 
<p V2 m p y/l V J 

N " = 44; (, - 2f) ~ 0(1) ' (50) 

N '"-T^M c "-" 2 + \ e )~ 0{cl) ' (51) 

where we have used the potential slow roll parameters 



ml fV'\ 

^^ 77 - (52) 



V 



v 

V'V" 



V = mp— > ( 53 ) 



( 54 ) 



Hence the non-linearity parameters for single field inflation (|45H46j) are given by 



Inl = !fa-2e), (55) 
t nl = (7]-2e) 2 , (56) 
gNL = || (2e ?? - 2T7 2 + e 2 ) • (57) 

Note however that we have not calculated the full bispectrum and trispectrum at leading order in slow roll, be- 
cause we assumed that the initial field fluctuations were Gaussian. If we included the contribution from the non- 
Gaussianity of the fields at Hubble exit, the bispectrum would have one extra term (|19p and the trispectrum would 



have two extra terms, (|36|) . The extra term for the trispectrum is at the same order in slow roll, because 12 1 
B As ^ki,k2,k 3 ) ~ 0(e3) and the second term of fl3BJ) is also of the same order. However Sccry, Lidsey and Sloth 
find [10], T ABCD (ki, k2, k/3, k^) ~ 0(1) which means the first term of (|36|) is suppressed by one less order in slow 
roll then the other three terms. However they find the contribution of this term is still too small to be observable, 
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even in the multiple field case [18[. All of these extra terms from the non-Gaussian field fluctuations are momentum 
dependent, while all of the non-linearity parameters are independent of momentum. 

In single field inflation, £ is conserved at all orders on superhorizon scales. Therefore no evolution of the bispectrum 
and trispectrum is possible after Hubble exit. Hence neither will be detectable in cosmic microwave background or 
large scale structure experiments, but from the 21 cm background there is the possibility that /jvl ~ 0.01 might be 
observable [27j |. 



2. Curvaton scenario 



In the curvaton scenario a weakly-coupled field (the curvaton field, x) which is light, but subdominant during 
inflation comes to contribute a significant fraction of the energy density of the universe sometime after inflation. After 
it eventually decays, it is the fluctuations in this field that produce the primordial curvature perturbation, £ j22l |. 

In general, the energy density of the curvaton is some function of the field value at Hubble-exit, p x cx g 2 (x*), an d 
hence the primordial curvature perturbation when the curvaton decays is of local form (|44p . In the sudden-decay 
approximation the non-linearity parameters are [Til . |28| 



Inl = 



4r 



1 



„/2 



5r 
~6~ 



(58) 



and |2i 



9nl 
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and tjvl satisfies (|48p. where the parameter r, is given by 
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,99" 



+ lOr + 3r' 



(59) 



3 Px 



3p x + 4p r 



decay 



(60) 



where p x is the density of the curvaton field and p r is the density of radiation and hence r satisfies < r < 1, . 

One can obtain significant non-Gaussianity if the curvaton does not dominate the total energy density of the 
Universe when it decays, r <C 1, in which case we have 



NL 



9nl 



T (1 + *? 

At \ g'- 
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(61) 
(62) 



One obtains a significant bispectrum, /«r > 1 for r <C 1 if gg"/g' 2 ^ — 1- On the other hand if gg" / g' 2 ~ —1 the 
bispectrum can be small even for r <C 1 [29j and the first signal of non-Gaussianity could come from the trispectrum 
through qktt. ~> 1 [23|, 



V. CONCLUSIONS 



In this paper we have given a general expression characterising the distribution of the primordial curvature pertur- 
bation due to an initial distribution of scalar field fluctuations during multiple-field inflation, using the cW-formalism. 

We have given expressions for the power spectrum, bispectrum and trispectrum including all terms at leading 
order in a perturbative expansion, allowing for scale dependence and cross-correlations of the fields. In particular 
the connected part of the primordial trispectrum consists of four terms each with a different momentum dependence. 
One term depends upon the intrinsic connected part of the trispectrum of the field fluctuations and one term depends 
upon the bispectrum of the fields. 

At lowest order in a slow-roll expansion the field fluctuations shortly after Hubble-exit are expected to be Gaussian 
and scale-invariant. In this case Lyth and Rodriguez [ll| showed that the bispectrum from multiple- field inflation 
can be parameterised by a single non-linearity parameter f^L dependent upon the second-derivatives of the local 
expansion, N, with respect to the field values, given in Eq. (|38p . The connected part of the trispectrum can be 
parameterised by two further non-linearity parameters, tjvl and gNL, given in Eqs. (T4"l"j) and (|4"2"|) . tjvl is another 
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function of the second derivatives of N [2l| . whereas gNL is a dependent on the third derivative. In the particular 
case where only one field generates the primordial perturbation we have tjvl oc ffj L [f| [24j|. However g^^ can give 
rise to a connected part of the trispectrum at the same order [23| and hence constraints on the primordial bispectrum 
do not necessarily constrain the primordial trispectrum in such models. 

In the case of standard single field inflation the trispectrum will be too small to ever be observed, but in alternative 
models such as DBI inflation [3(3| or the curvaton scenario [23| the trispectrum may be observable and could be an 
important test of such models. 

Note added: While writing up this work, similar results appeared on the arXiv [3lj . 
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